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ABSTBACT 

Exact,  explicit  expressions  are  obtained  for  the  field  of  a 
radiating  periodic  point  dipole  located  in  medltun  1  at  the  center  of 
a  spherical  shell  of  mediujn  2,  which  is  "bounded  "by  medium  3-  Tii© 
vedues  of  Z  ,  JX  ,  G    and  the  shell  radii  are  arbitrary.  The  reflected 
and  transmitted  fields  axe  examined  in  various  special  cases.   When 
the  shell  radii  become  infinite,  the  fields  agree  with  those  obtained 
by  fi.  K.  Luneberg  for  the  case  of  plane  waves  normally  incident  on  a 
flat  plate  of  medium  2  betv.'een  media  1  smd  3«  When  the  radii  are  in- 
finite and  media  2  and  3  identical,  the  Fresnel  formulae  for  normal 
incidence  of  plane  waves  on  a  plane  interface  between  two  half-infinite 
media  are  obtained.  When  the  radii  are  finite,  media  2  and  3  identical 
and  /A.^   =  ^        the  results  of  C.  T.  lai  are  obtained.  For  finite  radii, 

media  1  and  3  identical,  and  shell  thickness  small  compared  to  a  wave- 
length, the  fields  check  the  approximate  results  of  J.  B.  Keller  for 
the  fields  reflected  and  transmitted  by  a  thin  shell  of  any  shape, 
when  his  results  are  specialized  to  the  present  case. 


1. 


Introduction 

The  electromagnetic  field  due  to  a  radiating  periodic  point 
dipole   in  a  homogeneous  isotropic  medium  of  infinite  extent  is  well 
knovra  .     In  this  paper  we  intend  to  find  the  modification  of  the  field 
when  the  source  is  at   the  center  of  a  medium  "bounded  hy  a  spherical   shell 
of  different  material,  with  a  third  medium  of  infinite  extent   external 
to  the  shell.     It  is   expected  that  part  of  the  field  will  he  reflected 
from  the  shell  and  part   transmitted.     These  reflected  and  trsinsmitted 
fields   £ire  of  particular  interest.     This  prohlem  is  solved  exactly  and 
the  solution  is   examined  in  special  cases,   in  some  of  which  it  reduces 
to  known  results.     In  addition,    the  solution  is  used  to  test  the  accuracy 

of  a  recent  approximate  calculation  of  the  reflection  and  transmission  hy 

2 
an  arbitrarily  curved  thin  shell   . 

Formulation  of  the  Prohlem 

A  point   dipole  with  axis   in  the  z-direction  is  at   the  origin 
of  coordinates  in  medium  1,  with  constants   U-, ,   C'  t  »  and  O",  ,   and  hounded 

at  r  =  r     "by  medium  2.     This  medium  extends  from  r  =  r     to  r  =  r,    and  has 
a     '  a  0 

the  constants  ix.        £_,   and    (T-.     Bounding  the   shell  at  r  =  r.    is  medium 

3  of  infinite  extent   and  with  constants    U-,     fi_,   and   (T-,. 

The  electromagnetic   field  vectors  E  and  H  must  satisfy 
Maxwell's  equations   and  appropriate  boundary  conditions  at  the  two  inter- 
faces between  the  shell  and  the   surrounding  media.     In  addition  the  field 
must  vanish  at  infinity  and  also  contain  the  incident  field.      We  assume 
the  time  dependence  of  the  field  components  to  he  of  the  form  e         ,  where 
<*>  is  the  angular  frequency  of  the  dipole.     With  this  time  factor,   no 
volume  charge  present,    and  yu..    £f    aJid    <S'    each  constant,  Maxwell's  equa- 


1,    Stratton,  J.  A.;    Electromagnetic  Theory.   Mc&raw-Hill,   I9I+I,   New  York 
PP  U36-7.                                                    — 
Ramo,   S.   and  J.   R.  Whinnery,   Fields  and  Waves   in  Modern  Radio     Wiley 
19U4,  Hew  York,  p  14.30.         ~  ~" 

2.  Keller,  Joseph  B.  ,    "Reflection  and  Transmission  of  Electromagnetic 

Waves  by  Thin  Curved  Shells",   Washington  Square  College  Mathematics 
Research  aroup.  Research  Report  No.   I72-5,   Jan,19U«,  Report  prepared 
and  written  under  Contract  No,  W2g-099-ac-172  with  Army  Air  Forces. 


tions  "beooiae   : 


2. 


y.  E  =  0 

7.  H  =  0 

VxE  =  -^    H  (1) 


7x  H  = 

viae  re 


c 


e  -  £  —  (2) 


and  the  constants  have  the  appropriate  values  in  each  medium. 

The  "boundary  conditions,   requiring  continuity  of  the  tan- 
gential comoonents   of  S  and  H  and  the  proper  jtunp-discontinaity  of  the 

li 
normal   conraonents  of   the  fields,    are  given  "by  : 


at  r  =  r 
a 

\-^^» 

at  r  = 

^.  =  ^2. 

\  '  =3. 

'i^.  =  ^'A„ 

«2  ^^.  '  ^>3. 

^\. = /y-sn 

/^  =2„  ■=  /'3"3n 

(3) 


In  these  equations  E,      rex)resents   the  tangential  component  of  E  in  medium 

■'-t 

1,  E-  the  tangential  component  of  E  in  medium  2,  E,   the  normal  component 
'-■^  -'•II 

of  E  in  mediura  1,  etc. 
Use  of  Hertz  Vector 

We  introduce  the  Hertz  vector  A  which  satisfies 

V\..  ^±f^      A  =  0  (U) 

c 

It  is  easily  verified  that  the  field  vectors  will  satisfy  Maxwell's  equa- 

tions  if  we  define  them  in  terms  of  the  Hertz  vector  as 

3*  Ramo,  S.  and  J.  R.  Whinnery,  op.  cit.  p  I5U. 

'+.  Stratton,  J.  A.,  op.  cit.  pp  kSj)-h. 

5«  Ramo,  S.  and  J.  R.  Whinnery,  op.  cit.  p  I65. 


H  =    V  X  A 

(5) 

^-     I^    V(V.A)-     ^    A 

For  a  point  dipole  with  axis  in  th.e  z-direction,    in  a  medium  of  infinite 

extent,    the  solution  of  equation  (U)  is  fcund  to  "be  a  vector  with  a 

single  component 

-ikr 

A     =  A       _® (6) 

z         o  r  ^   ' 

where  A     is  a  constant  proportional  to  the  dipole  strength,   and 


k=  -f-v??  (7) 


Solution  of  the  Problem 


In  the  present  case  we  assume  the  solution  of  equation  (U) 
to  he  a  vector  with  the  single  non-zero  component: 


-Ikj^r 

ik^r 

\=K\ 

.A,!-^ 

0  <  r  <  r^ 

-ik^r 

ikgr 

A     -  A     ® 

*-3^ 

^a  <  ^  <  ^ 

z      ^Z       r 

-ik  r 

A     =  A     ^ 

r      <   T 

z       \       r 

r^  <  r 

(8) 


In  this   solution  the  constant  A     is  presumed  to  he  knovm  since   it  is 

determined  hy  the  dipole  strength,   and  tiae  constants  A,  ,  Ap,  A,,  A^  are 

to  he  determined.     Trom  eqiiations   (8)   it  can  he   seen  that  the  assumed 

solution  contains  the  incident  field  (6)  and  vanishes  at  infinity.     By 

suhstitution  one  finds  tliat   (8)   is  actually  a  solution  of  equation  (U), 
and  therefore  the  field  components,   which  are  ohtained  hy  using  (8)   In 
equations  (5).  will   satisfy  Maxwell's  equations   in  each  region.     These 


6.      ihid.  p  1+30 


u. 


field  components  axe: 
-ik,r 


^f.  ={  ^^"  ^'' [ '\  *  ?]^  ^^  ''  h^  *  7]]  F  ^^-® 


,     =[  a/^^2^  [ik^  .  ^]^  A3e'^2'  [.ik^  ^\]]l     Sine 


-ik^r 


[^,     =^Ai^e        3    ["ik^+ljll     sine 


COS  ' 


2c  --,8 


i  w  e'^r^ 


cos' 


i«*>e,'   r^ 


£ 


«. 


r      -ik^^r  r    ,     k  -,  ik,rr  k,        .-. 


(9) 


(10) 


, 1 smji 


£c 


£c 


I-ik  r  r  k  ,  ik  r   r     ^       k^  11 

|v  "''  i'-^,  ^  ^  -  7]  }  life  -- 


sinO  (l]) 


Hg.    =  Hr.     =  Evi    =0,  (i  =  1,2,3)  (l2) 

Frora  equations  (9)  and  (10)  we  see  that,  in  each  region.  H  .  and  E  are 
related  ty 

Hy  =  6'  E^(~^)   tmie  (13) 


5. 


We  must  now  determine  A.    (i=l,2,3i^)   so  that  the  field  com- 
ponents  (9)-(l2)  will   satisfy  the   appropriate  TDoundary  conditions   (3)at 
the  two  interfaces.      On  acco\int   of  equations    (12)  we  have  E^  =  Eft     , 
H^  =  H^  ,   and  of  course  E„  =  Er  and  H«  =  Hr   .     Prom  (12)  it  immediately 
follows   that  the  condition  on  the  normal  component  of  H  is  satisfied 
identically  at  "both  boundaries.      If  Hy    satisfies   the  condition  presented 
"by  equations   (3)  set  the  interfaces  then,  by  virtue  of  equation   (I3),  E 

will   satisfy  the  condition  for  En     as  given  in  (3).     Thus  the  field  com- 
ponents  (9)-(l2)  will  satisfy  Maxwell's  equations   and.  the  boundary  con- 
ditions provided  A.    (i»l,2,3,lf)  are  determined  by: 


oa     0  la  1  2a     2         3^     3 


0  = 


;    at  r=r^,  H^,  =  H^^ 
^2b  ^2  *  ^3b^3  -  \-b\     •   ^*  '^'b-  ^^-=  ^^* 


oa    0  la  i         2a     2        3^     3 


;    at  r=rg^,  E^^  =  Eg^ 


(lU) 


0  = 


°2b  ^2  "*■  °3b  ^3  -  ^-bh  '  ^*  """^b'  ^e."  ^«» 


v;here  the  follov;ing  abbreviations   have  been  used: 

-ik,r 


-ik^i 


B       = 
oa 

e 

\a  = 

ik^r 
1  a 
e 

"2a  = 

e-"='^ 

v  = 

e 

^2b  = 

e 

B 


3t> 


=  «il=2^b 


"la 

C       =  ^ 
oa             ^i 

'1 

r       =  e     -"■  ^ 

-^Va 

'^4  *}-\] 

L               a        r^J 

a         r 

&■' 

2       r^         ^2 
a  J 

ik2    -    ^2  _     i 

■2a               .    ■- 
^2 

^^2^a 

*^2 

~^^2'b 
0       -  e,, 

11.3  *  _  - . 

2       kg         i   ' 

^?] 

(15) 


\13 


3  ^ 


h  •  k] 


-^3^b 


c       =  e 

wii-u     ;7 


%-b 


^3 


6. 
(15) 


The  four  linear  equations  {Ih)   in  the  four  unknowns  A_ ,  A^, 
A   Aj^  have  the  solution: 

A 


(^l)St  -  ^3^''h-t^^^oep2.   -  ^2a°oa^  *  ^^b°2b  "  \^'h^^3.'o^-^c.^] 


V  -^  [(^t°3b  -  ^3t^l)^^Woa  -  ^oa°la)} 


(16) 


Ai, 


^  |(B3tC2t  -  B2^C3^)(B,^0i^  -  Bi^  C^^)| 


where 


^  =  [\^^3^-^3^^kt^^^lap2sr^2^ha^   *  (^^°2h-^2tC^l3)(23a°la-^laC3a)} 

The  values  of  A.  can  he  found  explicitly  "by  inserting  in  (l6)  the  ex- 
pressions given  in  (I5).  Thus  an  exact  solution  for  the  Hertz  vector  is 
given  hy  outting  into  equations  (8)  the  values  of  A.  given  "by  (I6),  and 

the  fields  are  given  "by  putting  these  values  of  A.  into  equations  (9)- 

(12). 

Discussion  of  the  Solution 

From  equations   (S),    (I5),    and  (lb)  v?e   ohtain  the  Hertz 
vectors   of  the  reflected  and  transmitted  fields.      They  are 


V  -4^1^a  *  H^^^a-4 


lk,r 


A     e         _       ,     e 
^  -r ^0  - 


-ik  r 


ilc^h  -ik_h 

e*^       F.M  +  e  J.P 

Le.  F.&  +  e  J.LJ 


(17) 


^u  — \  ^2  S  ~ 


-i[Va'*-  V^-'t^^ 


<l  .  s 

« 

ik  h               -ik_h 

e             F.G+e             J« 

L 

(18) 


v/here  the  new  symTaols   are  defined  "by: 


=HtvyN-^-y-'^K*yH*^a 


a  a  ' 


D         a 


The  Hertz  vectors  of  the  reflected  and  transmitted  fields  given  "by  equa- 
tions (17)  add  (18),  assume  simpler  forms  in  specisl  cases,  some  of  v;hich 
will  now  "be  examined. 


8. 


I,   When  li=0  or  r  =  r.    ttie  shell  of  medium  2   is  ahseut.   In 
a       D 

this  case  equations   (IJ)  and  (18)   'become: 


h 


ik^r 


-4Va  *  ^l^^a-^)] 


=  A, 


/       _/ 


Ir  ,  ,  (k,-k,)  ,         ,-1         r(k  -k   )      ,  ,  (£7-4)   1 


(20) 


-ik  r 
A     e ^  ^        e 


-i[k,r^+k3(r-rj] 


-2£ 


'd. 


^a  ^  ^a  J 

(21) 

If,  in  addition,  r  ■becomes  infinite  vdaile  r-r  remains  finite  and 
s  a. 

^-i^l^a 
•^0  r approaches  a  limit  A.^,    the  problem  reduces  to  that  of  a  plane 

wave  normally  incident  on  a  plane  interface  "between  two  half-infinite 
media.   In  this  case,  writing  r-r  =  x,  equations  (20)  and  (a)  redace  to: 


h'— 


A 


1^  —7 


(22) 
(23) 


9. 


Equations   (22)   and  (23)  axe  the  well-knoip  Fresnel  formulae  for   normal 
incidence^  and  equations   (20)  and  (21 )  are  generalizations  of  these 
fomralae.   If/*i  =/*3  equations   (20)  and  (21)  reduce  to  those  obtained  "by 
C.T.Tai     for  this  special  case, 

II,   When  r.    and  r     become  infinite  vdiile  r.-r_  =  h  and 
^  a     -i^l^a 

r-r     =  X  remain       finite  and  A     - — approaches  a  limit  A'    ,    the 

a  or  o 

problem  reduces   to  that  of  plane  waves  normally  incident  on  a  flat  plate 

separating  two   different  half-infinite  media.      In  this  case  equations 

(17)  and  (18)  become; 


ik,r 
A^    5 =  A'   e 

1.  T  0 


ik-x 


e^gC  ^1^3  *  ^3^1^  °°^  ^2^^  +  1(^1  ^'fl  +  ^2   ^1^3^  ^^  ^2^ 


(2U) 


-^V  -ik_(x-h) 

Av  ^ =  A«   e       ^ 

4         r  0 


^2  k2(  ^1^3  *  ^3^^   ^°^  ^2^  ■*■  ^^  ^i  ^3^2  ■*■  ^2\^3^   ^^^  ^2"^ 

Equations  (2l+)  and  (25)  are  the  same  as  those  obtained  by  R.  K.  Luneberg 

9 
for  this  ca.se'^, 

?•  Stratton,  J.A, ,  op.  cit.  p  1+92 

^»  Tai,  C.T.,  "Radiation  of  a  Hertzian  Dipole  Immersed  in  a  Dissipative 
Medium",  Harvard  University,  Cruft  Laboratory,  Research  Report  No,  21. 

°»  Luneberg,  R.X, ,  "The  Propagation  of  Electromagnetic  Plane  Waves  in 
Plane  Parallel  Layers",  Washington  Square  College  Mathematics  Research 
Group,  Research  Report  No.  I72-3,  June  19^7,  pp  U3-9.  Report  prepared 
and  written  under  Contract  No.  W28-099-ac-172  with  Army  Air  Forces. 


10. 


III.     If  h/r     «  1  and  |  kji     «  1  then  the  Herts  vectors 

given  hy  equations   (Ij)  and  (IS)   are  accurately  represented  by  the  first 
few  terms   in  their  Taylor's  expansions   in  powers  of  h.     The  first  terms 
of  these  expansions   are  given  helow  for  the  case  in  which  medium  3  is  the 


same  as  medium  1,  and     l/k,r    I  <<  1, 


ikj^r 


-i[2kir^  -  kir] 


=  A, 


^2^1 


(25) 


-ik,r 


,k^: 


-ik,r 


h' 


A 


k  ~r- 


=  A 


0       r 


^1"  "  ^~  liK  ^     f  i^2 


(26) 


By  using   (25)  and  (26)  in  equations    (9)-(l2)  we  obtain  the  reflected  and 
transmitted  field  components  expanded  in  powers   of  h.     These  are,  vjhen 


Tk^«i. 


ik,r 


H, 


-ik^r 
%    =  Ai^^-^r-(ikisine) 


^h^     ick^ 


V^ 


(-— i  sin  e) 


^  W£i 


E,    =  A^  ^-^r-(  -;^cose) 


-ikir 


e 


Wf^'r 


2c 


h 


^r  '  h  ~T^  ^tt:^  """^  ^^ 


oe'r 


^  =  ^e=\=^ 


(27) 

(28) 
(29) 

(30) 

(31) 
(32) 

(33) 


11. 


In  these   equations  the  letters  H  and  T   in  the  subscripts   denote   reflected 
and  transmitted   field  components,    respectively. 

Comparison  vith  an  Approximate  Solution 

An  approximate  solution  of  the   reflection  and  transmissioii 
of  plsme   or  spherical    electromagnetic  waves   from  an  arbitrarily  carved 
thin  shell  has   recently  "been  given  "by  J.  B.   Keller      .     The  results,  when 
specialized  to    the  case  of  plane  waves   incident   on  a  thin  flat  plate, 
e^ree  with  those  of  R.   K.    luneberg     .      In  order   to  obtain  another  check 
on  the  approximate  solution  we  will   specialize  it   to  the  case  of  a  point 
source  at   tlie  center   of  a  spherical    shell.      The  fields  obtained  in  this 
way  should  agree  v;ith  those  given  by  equations   (27)-(33)   above. 

When  equation  (U5)   in  the  paper  by  J.   3.   Keller  and  the 

corresponding  equation  for  H  are  specialized  to  this  case  it   is    found 

that  the  reflected  and  transmitted  H  ,   B.q  ,    and  E^are  all  zero.     We  also 

find  that  the  reflected  and  transmitted  H^    and  Eq    axe  given  by  equations 

(27)-(3'^)«     Il^e  reflected  and  transmitted  E     components  become: 

'^1"  2ck, 

Er     =  A,  «__  (   -  ^  cos  e)  ^  (34) 


'^S 


\ 


A^H.(A^.AJI- 


0     r 
a 


e     ^      ,2ck^ 


2ci 
^ZTT^  ^°^  ^  )  (35) 


'UJZ 


f 


We  thus  see  that  the  approximate  solution  agrees  with  the 
results  of  the  special  form  of  the  exact  solution  given  by  equations  (27)- 
(33)  in  all  components  except  E^  and  E^  .   Since  equations  (27)-(32)  are 
valid  only  i^hen  \  k-^r^  |  »  1  and  |  k^^r  [  >>  1,  the  ^proximate  solution  can 
give  good  results  only  \inder  the  same  conditions.   The  first  condition  was 
to  be  expected  from  the  derivation  of  the  approximate  solution,  which 
assumed  that  the  source  v/as  many  wave-lengths  from  the  obstacle.  In  that 
determination  it  was  also  assiimed  that  there  v/as  no  singularity  in  the  re- 
flected  field.   In  the  special  case  here  considered  there  is  a  singularity 

10.  See  reference  in  footnote  2  on  p  1, 

11.  See  reference  in  footnote  9  on  p  9. 


12. 


at  r  =  0  and  hence  it  is  to  "be  expected  that  the  result  will  hold  only  at 
a  distance  of  several  wave-lengths  from  the  singularity,  i.e.,  1  k^r  )  »  1, 

In  the  case  of  the  components  E   and  E   ,  equations  (3^)  and 
(35)  hecowe  identical  with  equations  (3I)  and  (32)  when  r  =  r   hut  differ 

for  other  values  of  r.  The  error  in  the  radial  dependence  of  E   is  due 

to  the  sensitive  dependence  of  E  on  the  method  of  evaluating  an  integral 

occurring  in  the  approximate  solution. 

From  the  ahove  comparison  of  the  approximate  and  ezact  solu- 
tions it  can  he  seen  that  the  approximate  solution  is  correct  except  for 
the  longitudinal  field  component  E^,  and  except  ne?r  the  singularity  of 

the  reflected  field.   Since  the  occurrence  of  singularities  in  reflected 
fields  is  uncommon,  as  is  also  the  occurrence  of  longitudinal  components, 
the  comparison  indicates  that  the  approximate  solution  should  usually  he 
acceptable. 
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